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Abstract 

We study the notion of stochastic stability with respect to diffusive perturbations for flows 
with smooth invariant measures. We investigate the question fully for non-singular flows on 
the circle. We also show that volume-preserving flows are stochastically stable with respect to 
perturbations that are associated with homogeneous diffusions. 
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1 Introduction 

The notion of stochastic stability of the invariant measures of a continuous dynamical system 
was formulated as early as 1933 by Pontryagin, Andronov and Vitt [T7j and Bernstein [3J. They 
solved the problem completely for flows on the real line that possess a globally attracting set of 
equilibrium points. Such flows of course are the gradients of some potential function. Freidlin 
and Wentzell [11] extended the classical work to M. n , again with globally attracting equilibrium 
points, by defining dynamically a quasi-potential associated with the flow. In the case of a unique 
equilibrium, stochastic stability is proved and a formula for the perturbed invariant measure was 
produced in terms of the quasi-potential. Stochastic stability for flow on a compact manifold was 
discussed in [14J, |16| and extended to the case that the maximal invariant set of the flow consists 
of a finite number of hyperbolic basic sets. Since then, it seems that very little attention has been 
given to this problem. 

In the meantime, stochastic stability of invariant measures has been a topic of continuing interest 
and surprising difficulty for discrete time dynamical systems theory. The main directions where 
progress has been made are: structurally stable systems [15], hyperbolic systems [13], [15] . [18] and 
unimodal maps [3], [2]. 
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In this note we study the question of stochastic stability for the (overlooked) case of non-singular 
flow on the circle and for a volume preserving flow on a manifold. 

We will let h : M — > TM be a smooth vector field on a compact Riemannian manifold M with the 
associated differential equation in local coordinates: 



By standard results, h generates a smooth, global solution flow $':lx¥-> M. All equations in 
this note will be assumed to be written in local coordinates. 

We will denote by d the distance on M, by B the Borel cr-algebra on M, by m the normalized 
Riemannian volume (Lebesgue measure) on M, and by M the space of Borel probability measures 
on M with the weak topology. We will use — to denote weak convergence in Ai. 

Consider a vector field h for which there is a unique physical, ergodic measure fiQ G M,, An 
ergodic measure /j,q is physical if for m-a.e. x G M, u X) t '■= y f S$t( x ) Ho a s n — > oo, i.e. for all 
(/) G C°(M) we have £ jj <j>(&{x)) -»• / (pd^o,. We will also assume that no is absolutely continuous 
with respect to the Lebesgue measure on the manifold and has a smooth density function po. 

By stochastic stability we mean stability of no under small stochastic perturbations of (/>*, which we 
now define. Consider the stochastic differential equation 



which we interpret as a Stratanovich integral equation on the tangent space of the manifold M. 
Equation fll .2 j) is associated with a diffusion process on M. For general background on diffusions 
on a manifold see for instance |12|. Chp. V], [8] or [9]. We call Q1.2|) a small stochastic perturbation 
of (II. ip . A class of perturbations refers to a collection of equations of the form (|1.2j) with 7(x) 
chosen from some class of nonsingular matrix fields. By stochastic stability with respect to a class 
of small stochastic perturbations, we mean that each flow defined by (II. 2ft within the class has a 
unique ergodic stationary measure \x e and /i £ ^ /io as e -> 0. 

Note that any weak limit point [i of the set {/u e } as e — )• is called a zero-noise limit measure for 
h. It is well-known that all zero- noise limit measures are invariant under h [2]. See also [6] for a 
recent related work using zero- noise limits. 

We note here that no counterexamples to stochastically stability for physical ergodic measures of 
dynamical systems, except one in which the "random" perturbations have a distinctly non-random 
character [I]. 

By standard diffusion we mean for each x, ~f(x) is a multiple of the identity matrix. 

Also note that we could consider larger classes of perturbations that include both stochastic and 
deterministic perturbations by adding a term eg dt to the right hand side of (|1.2p , where g is selected 
from some class of vector fields. 





dx = h(x) dt + ->/ej(x) o dW, 
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We can associate with (jl.ip and (|1.2p a generator 

£ e = fc.V + |r:W (1.3) 

where T(z) = r ){x)^{x) T and ":" denotes the double inner product of matrices (the inner product 
consistent with the Frobenius norm). The formal L 2 adjoint of C is 

C* = -V-/i+|V- V-T. (1.4) 

Note that T(x) is a symmetric positive definite matrix field. For standard diffusion, T(x) is the 
identity matrix and zC e and C* reduce to: 

C e = h ■ V + |A and C* t = -V • h + |a. (1.5) 

Here A is the Laplace-Beltrami operator on M. 

Probability density functions evolve under (|1.2j) via the Fokker-Planck equation 

I = < L6 > 

(For e = this is usually called the Liouville equation.) Suppose that po is an invariant measure 
for (jl.ip and the density p^ of /xq is a C 2 function. Then £gPo = 0. Further, if p e is the density of 
a stationary measure /i e for (jl.2p . then 

Oe = 0- (1.7) 

Since M is compact and T is positive definite, the density p e for the perturbed system (|1.2p exists 
[7] and it is unique, nonzero and smooth [12, Prop. 5.4.5]. 

Now we wish to treat the behaviour of p e as a perturbation problem. Define r e (x) = p e (x) — po(x) 
so that 

p t = Po + r e - (1-8) 

It follows that r t exists, is unique and is smooth because p t exists and is unique and smooth, and, 
po is unique and smooth by assumption. Since J p e = 1, then any solution r e must satisfy: J r e = 0. 

Substitution into (jl.7p gives a perturbed elliptic equation with a constraint: 

e - V • V • (Tr 6 ) - V • (hr t ) = ~ V • V ■ (Tpo), subject to ^ r e = 0. (1.9) 
In the case of standard diffusion, this simplifies to: 

| Ar e = V • (/ir e ) - | Ap , subject to J r e = 0. (1.10) 
To complete this approach we need only show that r e goes to zero as e goes to zero. 
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Note: 

Instead of (jl.8p we could suppose that p e has the form 

p e {x) = po(x) + eu{x) + r e (x) (1.11) 

for some smooth function u(x) and remainder r e (x) of order o(e). However, this is a stronger 
assumption than we need and may impose unnecessary constraints. In one dimension one can solve 
directly for u, but in higher dimensions a solution seems unlikely. 

2 Two special cases 

In this section we point out that stochastic stability with respect to standard diffusion follows easily 
from known results in two extreme cases, volume preserving flows and gradient flows. 

First suppose that h is volume preserving and the perturbation is by standard diffusion. According 
to \12\ Theorem 5.4.6] the measure cm (c constant) is an invariant measure of a diffusion generated 
by A + h if and only if h is volume preserving. Rescaling (jl.2p by 1/e implies that the volume itself 
is the invariant measure of the perturbed process for any e > 0. This gives us the following: 

Proposition 2.1 If h is volume-preserving and the volume, m, is ergodic, then m is stochastically 
stable under standard diffusion. 

We will extend this result slightly to the case of homogeneous diffusions in § 4. 

Secondly, consider the case that h is a gradient flow, i.e., there exists a smooth, real-valued function 
H on M such that h = — gradi?. This case is the direct generalization of the early results [17], [4] 
and a special case of result in [11]. In \\.2\ Theorem 5.4.6] we find that if h = — gradff then the 
diffusion generated by A + h has as its invariant measures multiples of exp(— 2H{x)) dx. Rescaling 
by 1/e implies that the invariant measure for 11.21 is: 

c e e~ 2H{x)/e dx. (2.1) 

This measure converges as e — > to a measure supported on the set of minimum points of the 
function H(x). In the simplest case we have: 

Proposition 2.2 If h is a gradient flow with potential H , and H has a unique minimum point, 
then the delta measure on this point is stochastically stable. 

If there is not a unique minimum, then the delta measure at each of the minimum points is ergodic 
and any linear combination of them is (nonergodic) invariant. Freidlin has studied the asymptotics 
of transitions between the local minima as e — > [10] , 
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3 Flow on a Circle 

Suppose M = T 1 and h is smooth and nonzero. Then h generates a smooth flow and possess a 
smooth ergodic measure supported on the entire circle, given by the density po(x) = c/h{x). 

On the circle, (11. 9h becomes: 

| (Tr e )" - (hr e )' = -| (Tp )" subject to J r e = 0. (3.1) 

For e = 0, note that there is a one dimensional space of solutions of (13. lft : ro = c/h{x) is a solution 
for any constant c. If we restrict to J ro = 0, then ro = is the unique solution. As stated before, 
the existence, uniqueness and smoothness of the solution r t of (I3.ip follow from the assumptions 
on po and the existence, uniqueness and smoothness of p e . However, there is also an elementary 
proof of these facts in this case. 

Assume that h, po and T are all at least C 2 smooth on the circle. Then r e is a solution of (j3. 1 j) if 
and only if it is a solution of 

e -(Tr e y-hr € = - e -(r Po y + C e , (3.2) 

for some constant C £ . 



Proposition 3.1 If r e is a solution of i3.1\) . then r e converges to zero in the L 2 norm as e — > 0. 
Proof: 

We will assume without loss of generality that h(x) > in local coordinates. Let r e be a solution 
of (|3.ip subject to f r = 0. Multiplying equation (|3.2p by r e and integrating over T 1 we have: 

Here and subsequently, J will mean integration over the circle, T . Note first that the final integral 
f C e r e dx is zero by constraint. Rearranging and successively integrating by parts we obtain: 



hr 2 dx 
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i J Yr t r' e dx + ^ J (Tpo)'r e dx 

iI T i r - dx+ U iTnYr - dx 

J T'r 2 dx + ^ j (rpo)V 6 dx. 



Thus we have 



J(h- e l T')r 2 dx= e - J(Tp o yr e dx. 
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Denote a = min h. In the case that T is not identically constant, note that we have max(r') > 
and so 

h-~V>a/2 (3.3) 
provided that e < 2a/ max(r'). (If T is identically constant, then (|3.3p holds for any e.) Therefore, 

| J r 2 e dx<^ J \(Tp )'r e \dx. 

Applying the Schwartz inequality to the last integral we have 

flNII< IlKrpoyiblKlh. 

or 

||r e || 2 < -\\(Tpo)%. (3.4) 
a 

Thus we have in fact that ||r e ||2 = 0(e). □ 

Next we show that r' e = 0(e) in the I? norm. Multiplying equation (|3.ip by r' e and integrating 
over T 1 we have: 

\ J{Tr e )"r' e dx- J{hr t )'r' e dx = - e - J (Tp Q )"r' e dx. 
Expanding the derivatives and integrating by parts we obtain: 

J(h- |r')(r0 2 dx = -J(ti- e -T")rJ e dx+ € - J (Tp )"r' e dx. 

Note that for e < 2a/(3maxr') (or any e in the case T' = 0) we have h — > a/2, and so 

J(r' e ) 2 dx< J \h> - e -T"\\r t r' e \dx+ e - J \(T Po )"r' e \ dx. 



o 
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Let 

P = m^{\h'\ + - -|r"|). 

3 max I ' 

(In the case V = 0, we may take j3 = max(Zi'). Then using the Schwartz inequality, we obtain: 

lK||i<^||r e || 2 ||^|| 2 + l||(rp )"||2||<|| 2 . (3.5) 
a a 



Using (|3.4|) we have: 

1 1 r 11^ < a 



^ll2<ef^ll(rpo) / l|2 + -||(r/ 3 o)"|| 2 

or a 



Thus we have: 

Proposition 3.2 As e goes to zero, \\r' t \\2 is of order epsilon. 
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Now, by the Poincare-Wirtinger inequality in one dimension [51 p. 146.], the L°° norm of r t is also 
of order 0(e). (An elementary proof of this also exists in this case.) Since r t is smooth and the 
circle is compact, the following proposition holds. 

Proposition 3.3 The solution r e of h3.1\) converges to zero uniformly as e — > 0. 
In terms of Stochastic Stability we have shown the following result. 

Theorem 3.4 Let h € C 1 be a nonsingular flow on the circle with an absolutely continuous in- 
variant measure fio with density po € C 2 . Then fio is stochastically stable with respect to any 
perturbations in the class C 2 . 

4 Volume-Preserving Flows 

Next consider smooth vector fields h that preserve volume, thus V • h = 0. If volume is also 
ergodic for the flow defined by h, then the constant multiples of the volume are the only absolutely 
continuous invariant measures. 

In this case (|1.9|) becomes: 

€ EC f 

- V- V • (IY e ) - V • (hr e ) = -— V • V • V, subject to / r £ = 0. (4.1) 
We note that V • h = implies: 

V • (hr e ) = h-Vr t 

If T is a multiple of the identity, or, if T is constant with respect to x (homogeneous), then (|4.ip 
simplifies to: 

V • (hr e ) = | Ar e , subject to J r e = 0. (4.2) 



Proposition 4.1 The only smooth solution of \4-<ty is r e = 0. 
Proof: Multiplying both sides of the equation by r e and integrating we have: 

J (h ■ Vr e )r e = | J r e Ar e 



J (hr e )-Vr € = ~\\ |Vr e 



(4.3) 
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On the other hand, r e is a smooth solution of (|4.2p if and only if there is a smooth divergence- free 
vector field v e such that 

hr e = - Vr e + v e . 
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We may take the inner product of both sides of the equation with Vr e and integrate to obtain: 




We can conclude from (|4.3I) and (I4.4p that both J (hr e ) • Vr e and J |Vr e | 2 are zero. This along with 
the constraint J r t = implies the result. □ 

We note here that there are other simple proofs that J (hr e ) ■ Vr e = 0, including a dynamical 
argument. To see this, note that (hr e ) ■ Vr e = (h ■ Vr e )r £ is in fact the directional derivative of 
r 2 along the vector field h. If we integrate this expression along a trajectory, ergodicity and the 
assumption that r t is bounded imply the claim. 

The implications for stochastic stability are the following: 

Theorem 4.1 If h preserves volume, then the volume is the unique zero-noise limit under any 
homogeneous diffusion. If volume is ergodic, then it is stochastically stable with respect to homoge- 
neous diffusion. 



5 Discussion of non-smooth invariant densities 

In general, it is rare for p,Q to have a smooth density, rather, it is usually supported on an attractor 
which has dimension less than that of the ambient space. Further complications include that the 
physical measure po will not be the unique ergodic measure supported on the attractor. 

If the method could be generalized it might be as follows. The Liouville equation must be considered 
as acting on distributions. The densities for the SDE (II. 2p however remain smooth. With these 
considerations, (jl.8p cannot be the right assumption. Rather, we need to show something of the 
form 

p t = A e p , 

where A e is a smoothing operator that acts on distributions and approaches the identity as e — > 0. 
It is natural to assume then is that A e is itself generated by a diffusion, i.e. 

A e = e £L * 

where L = T(x)A and T must be determined. In this circumstance, with standard diffusion as the 
class of perturbation, C*p € = becomes: 

V • h(A eP0 ) = -A{A ePo ), 
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and the problem is to find an appropriate f (x). 

Acknowledgment: The authors wish to thank Martin Hairer for pointing out important references 
and making several helpful comments. 
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